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Abstract 

Productions and decays of spin-singlet S, P— wave charmonium states, r\ c and h c ( x Pi), 
in the e + e~ annihilations are considered in the QCD multipole expansion with ne- 
glecting nonlocality in time coming from the color-octet intermediate states. Our 
approximation is opposite to the Kuang-Yan's model. The results are B(ij) — > 
h c + 7T°) « 0.3%, B(ip' -> ri + 7) w 0.34%, T(r]' c -> J/V> + 7) » 0.26 keV and 
T(h c -> J/^ + vr°) « 2.5 keV. 

PACS numbers: 14.40. Gx, 13.25. +m, 13.40.Hq, 11.50.Li 
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Among the heavy quarkonium states, the spin-singlet P— wave state (h C) bCPi)) is 
extremely difficult to be produced in the e + e~ annihilations, because it has J PC = 
Such a state in the charmonium family was first discovered in pp annihilation 
by E760 Collaboration at Fermi Lab in 1992 (almost twenty years after the discovery 
of charm quark) with the subsequent decay h c — > J/ipn° JTJ : 

T(h c ) ~ 800 keV, 
T(h c -> J/tp + 7T°) ~ 0(1) keV, (1) 
T(h c -> J/ip + tttt) < 0.18 T(h c -> J/ip + n ). 

The data supports Voloshin's approach || instead of Kuang-Yan's model || as dis- 
cussed in Refs. [0 and [|J]. Another spin-singlet S-wave state ?7 C can be reached in 
the e + e~ annhilation via a spin-flip radiative decay ip — > ?7 C + 7- However, this state 
has not been confirmed yet. Various potential models predict m(rj' c ) to be around 
« 3.6 GeV. 

Since BES (Beijing Spectrometer) is planning to search for h c and rf' c in the decay 
channels ip' — *> 77^7, h c ir° and h c Tnr, it is worth while to estimate the branching ratios of 
these decays in the frameworks of QED and QCD multipole expasions, and compare 
the results with predictions of the Kuang-Yan's model as well as with the data from 
BES and other experiments. This is precisely the purpose of this report. 

First, we discuss El and Ml raditaive transitions of charmonium states, ip — > 
Xcj + 7, Xcj J/ip + 7 and ip — > i] c + 7, in order to extract various matrix ele- 
ments such as (lP\r\nS) with n = 1,2 and (lS\r 2 \2S) from the experimental data. 
These matrix elements enter various hadronic transitions between charmonium states 
that will be considered later in this work. We also estimate the branching ratio for 
■0 ~ * Vc + 7- After discussing radiative transitions of charmonium states, we briefly 
recapituate how to describe hadronic transitions between heavy quarkonia in QCD 
multipole expansion and make a simple approximation which ignores nonlocality in 
time due to the color-octet intermediate state between subsequent emissions of glu- 
ons. One can estimate the size of the color-octet Green's function (G%(m^i)) from 
the measured decay rate for ip' — > J /ip + tctc (or rj). For this purpose, one has to use 
the matrix element (15'|r 2 |2S') determined from ip' — ► r/ c + 7. One can use this G§ to 
compute other decays such as h c — > J/tp + 7r°(or tttt) and ip' — > h c + n° (or tttt). Our 
results are summarized at the end. 

It should be emphasized that one cannot make precise predictions for decay rates 



1 



of hadronic transitions between heavy quarkonia because of our ignorance of color 
confinement in QCD. Any attempts to estimate such quantities are admittedly model 
dependent. Kuang-Yan's model takes into account the nonlocality of G$ in time, 
using the string excitations near the flavor threshold. Our approximation will be 
opposite to Kuang-Yan's model in the sense that we ignore nonlocality of Gg in time, 
and regard Gs as a constant. Predictions of both apporaches are a kind of rough 
estimates, and the actualty may lie somewhere between two approaches. In this 
sense, our predictions can be considered as complementary to Kuang-Yan's model. 

In order to describe hadronic transitions between heavy quarkonia, we have to 
know such matrix elements as (nS\r\lP) and (lS\r 2 \2S) . In previous estimations, 
these quantities were replaced by the radius of the initial or the final quarkonium. 
However, such estimates can be improved, since these matrix elements are also rel- 
evant to electric dipole and spin-flip radiative transitions bewteen heavy quarkonia. 
In this work, we work in the nonrelativistic quantum mechanics, neglecting the rela- 
tivistic corrections, LS coupling and hyperfine splittings. Therefore, the results may 
be affected by a factor of ~ 2 due to such neglected effects. 

Let us first consider electric dipole transitions, whose decay rates are described by 

T(E1) = A« e\ cv s S lf (2J, + 1) \(nS\r\mP)\ 2 . (2) 

Here, J/ is the spin of the final quarkonium, and to is the energy of the emitted 
photon. Sif = 3 for h^Pi) — ► r? c + 7 and r\ — > h c ( l Pi) + 7, whereas Sif — 1 for other 
transitions. 

Comparing (fj) and the measured decay rates of XcjO-P) ~^ J/ip + *y and ip' — > 
X c j(lP) +7, we get 

|(LS|r|LP)^| » 1.58 GeV^ 1 = 0.3 fm, (3) 
|(2S|r|LP)^| w 1.97 GeV" 1 = 0.4 fm, (4) 

where we have taken the weighted average over J = 0, 1, 2 states. Using m(h c ) = 3526 
MeV and taking m(r/ c ) = 3600 MeV tentative value, we get 

r(/ic -> Vc + l) ~ 460 keV, (5) 
r(r}' c ^h c + j) w 6.5 keV, (6) 

Since |(nS'|r|lP)^|'s with n = 1,2 are determined from the experimental data on 
T(Xcj — > J/ip + 7) and r^' — > x c j(1P) + 7), our predictions (H) and (6) are inde- 
pendent of specific potential models. 
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Next, we consider a magnetic dipole radiative transition, n 3 Si — > m 1 ^ + 7. The 
decay rate for such a spin-flip radiative decay is given by 

r(y/-tfc + 7 ) = \aQl— 2 \{lS\ 3 ^r/2)\2S)\ 2 (7) 

- 4 k^i^w- ( 8 ) 

3 24 

In (7), j'o(^) is the 0— th order spherical Bessel function, jo(x) = sinx/x. We have 
used the long-wavelength approximation (ur << 1) in (8). From the measured decay 
rate — > ?7 C + 7) = 0.78 ± 0.19 keV, one can extract 

\(lS\r 2 \2S}\ ^2.55 GeV~ 2 (9) 

which satisfies the bound (20) derived from (19). The main point of this work is the 
following : one can get the matrix element |(lS'|r 2 |2S')| from ip — > ^ c + 7. This matrix 
element will be used later in order to extract G% from ip — > J/ip + 77, as described in 
the following. 

Now, let's consider the possibilty of finding r\ c state via a radiative transition 
ip' — > 77^+7. Assuming m(^) = 3.60 GeV, the energy of the emitted photon is u; = 84 
MeV, and the matrix element in (7) can be approximated as 1. T(ip — > ^+7) ~ 0.94 
keV, or 3.4 x 10~~ 3 in the branching ratio. Thus, one would get ~ 340 rj^s among 
10 5 ip r, s. The subsequent decay of rj c is mainly hadronic through annihilation into two 
gluons : T(rj' c — > 2g) « 1 — 4 MeV. The decay rate of the spin-flip radiative transition 
r)' c — ► J/-0 + 7 is about 0.26 keV, using the matrix element | (lS\r 2 \2S) | which was 
determined above, (|]). 

Let us now discuss our main subject, the hadronic transitions between heavy 
quarkonium states in QCD multipole expansion @ @ @. The basic formulae can be 
found in previous works, and we write down the relevant equations only here. The 
decay rates for spin-flip transitions such as l x Pi — > l 3 Si + 7r° and n 3 Si — > l x Pi + 7r° 
are given by |J 

T(1 1 P 1 - l 3 ^ + vr ) = F(n 3 S 1 -> l 1 ^ + tt°) = J- (A J S p) 2 |&|. (10) 

Here, A = 1.7xl0~ 3 GeV 3 is one-third of the matrix element (vr°|7ra s E a (0)-B a (0)|0), 
and the matrix element I$p is defined as 

2V3 
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(nS \Gs(E)r + rG 8 (E)\lP), (11) 
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where Gg{E) is the Green's functions for the color octet QQ intermediate states : 

Gs(E) = £ (12) 

To get the absolute decay rate for hadronic transitions between quarkonia, it is im- 
perative to know more about "isp" defined in ([□]). However, it lies beyond our 
ability to calculate I$p from the first principle in QCD, because of our ignorance of 
the confinement in QCD. Therefore, we have to make some reasonable approxima- 
tions. We take an opposite limit to the Kuang-Yan's approach, namely that we ignore 
nonlocality in time coming from the color-octet Green's function, G$(E). We simply 
assume that Gg(E) is a constant with dimension of inverse mass as in Ref. ||], since 
this assumption does not lead to any contradictions to exsiting data : 

G$(E) = G$ = constant. (13) 

Under this assumption, we have 

4v/3 
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G 8 (nS\r\lP), (14) 



and one can use the informations on the matrix elements on r, Eqs. (3) and (4). 
Taking m c = 1.65 GeV and using (3) and (4), we get 

T(h c (lP) -> J/ip + 7T°) = 0.095 (—^—) |G 8 (GeV- 1 )| 2 keV, (15) 




r(V' - h c {lP) + tt°) = 0.032 ( ] |G 8 (GeV- 1 )| 2 keV. (16) 

\m c (GeV) / 

Note that (|T5|) is smaller than the measured value (1) by an order of magnitude if we 
assume G$ ~ 1 GeV^ 1 , as in Ref. ||. 

In order to get informations on G&(E), let us consider another hadronic transition, 
ip' — > J/ip + rj, which occur through El — El and El — M2 transitions in QCD 
multipole expansion [H] |5|. Its amplitude is proportional to the following matrix 
element 

Iss = \{mS\r l G{E)r l \nS) 1 (17) 
which is simplified to the following under the assumption (13) on Gg(E), 

Iss = Ig s (mS\r 2 \nS). (18) 
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Therefore, once (2S\r 2 \lS) is known, the absolute decay rates for these decays can be 
readily obtained. However, potential model calculation of this matrix element is not 
available on the contrary to dipole matrix elements (/ | r | i) and the mean square 



radius of a quarkonium, (i 
sum rule to derive 



i). In Ref. M, I have used a quantum mechanical 
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\(1S \ r 2 \ nS)\ 2 < — . 

m Q \^nS - &1S) 

From |(15'|r 2 |15')^| = 4 GeV~ 2 , the bound (19) was found to be 

|(l,S|r 2 |2S)^| 2 < 16.5 GeV" 4 . 



(19) 



(20) 



By considering if) — > J/ij) + r\ and using (20), one obtains a lower bound on G 8 || : 

|G 8 | 2 > 10.4 GeV~ 2 . (21) 

Assuming G%{m^/) ~ Gg(m hc ), we get the following lower bound on h c — > J/ip + n° 
from (15) : 

T(h c -> J/ij) + 7T°) > 1.0 keV, (22) 

This is consistent with E760 data ([I]), and moreover, very close to it. Of course, this 
would be uncertain by a factor of ~ 2, depending on the choice of m c . (In Ref. 0, the 
mass of h c was assumed to be 3.51 GeV, and thus numerical results obtained there 
are slight different from the results in the present work.) 

One of the main points of this work is that one can actually improve the lower 
bound (22) obtained in Ref. Q, using the estimate on (lS\r 2 \2S) obtained from 
■0 ~" * Vc + 7- From @ and — > J/ip + rj) = (7.51 ± 1.41) keV, one obtains 



26.4 GeV 



-2 



(23) 



Using this \Gg\ in (17) and (18), we finally get 

T(h c -> J/^ + 77°) ss 



2.5 keV, 
0.84 keV, 



(24) 
(25) 



Note that our prediction for B(ip' — > /i c +7r°) = 0.3% is very close to the current upper 
limit, and should be checked in the near future at BEPC. Also, these two predictions 
in our approach are remarkably in accord with predictions by the Kuang-Yan's model 
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with «Af ~ 10a e- However, we have only one gauge coupling (namely, oim = cue) in 
our model, and this agreement of two approaches must be regarded accidental. 

Another hadronic decay mode, l x Pi — > l 3 Si + nn, does not receive any contribu- 
tion from the trace of the energy-momentum tensor in QCD, and is not enhanced 
over 1 1 P 1 — > l 3 ^! + 7T° 0. (There is no anologous decay ip' — > h c + mr because the 
phase space is not avaiable.) Using the result of Ref. f2|, we predict 

T(h e ^J/1> + *°) " 30 " °' 16 ' (26) 
Here, A measures the gluonic contributions to the energy-momentum of a pion. This 
is consistent with E760 data, and very close to the current upper limit for A ~ 2. 
Therefore, h c — > J/ipim may be observed in the near future. 

A remark is in order. As mentioned before, we have made an approximation 
treating the color-octet Green's function G$(E) as a constant, and determined it from 
the measured decay rate for ip — > J/ip + rj. One can actually check the consistency 
of this approximation making the following observation. The ratio 

^r^? = 2A X 1(r * I^OV-lP (27) 
"(Ac -> r/ c + 7) 

depends on \G$\, which can be used as an independent determination of |G§| from the 
measurements. Another ratio between T(ip' — > /i c + 7r°) and r(^' — > x c j(1P) +7) can 
be useful as well. Thus, our assumption on Gg is simple, does lead to predictions on 
various hadronic transitions between heavy quarkonia, and can be checked through 
the measurement of the ratio (27). 

In conclusion, productions and decays of r\ c and h c in the e + e~ annihilations 
through i\j decays are considered in the framework of QCD multipole expansion, 
assuming the Green's function for the color octet states is a constant. Using the matrix 
element | (15'|r 2 |25') | extracted from the spin-flip radiative transition ip — > r) c + y, we 
could estimate the size of |G 8 | 2 by fitting ip — > J/ip + 77. This leads to predictions 
of absolute decay rates, B(ip' -»• t/ c + 7) « 0.34%, T(?7^ -> J/ip + 7) w 0.26 keV, 
and (24), (25), which improves the lower bounds on the decay rates for the same 
processes obtained in Ref. |§]. We also have suggested that our assumption ( |i~3|) 
could be checked by measuring the ratio (27). We hope high statistics experiments 
at BEPC and other places searching for rj c and h c can test predictions made in this 
work. Having ~ 10 5 ip decays would be enough to get the spin-singlet rf' c and h c 
states and their subsequent decay properties. 
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